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Voltage-modulated magnetism in magnetic/BiFe0 3 heterostructures can be driven by a combination of the 
intrinsic ferroelectric-antiferromagnetic coupling in BiFe0 3 and the antiferromagnetic-ferromagnetic 
exchange interaction across the heterointerface. However, ferroelectric BiFe0 3 film is also ferroelastic, thus 
it is possible to generate voltage-induced strain in BiFe0 3 that could be applied onto the magnetic layer 
across the heterointerface and modulate magnetism through magnetoelastic coupling. Here, we 
investigated, using phase-field simulations, the role of strain in voltage- controlled magnetism for these 
BiFe0 3 -based heterostructures. It is predicted, under certain condition, coexistence of strain and exchange 
interaction will result in a pure voltage-driven 180° magnetization reversal in BiFe0 3 -based 
heterostructures. 

■ t is accepted that voltage-modulated magnetism in magnetic/BiFe0 3 (BFO) layered heterostructures is based 
I on a combination of intrinsic coupling between the coexisted ferroelectric and antiferromagnetic orders in 
I BFO, and the antiferromagnetic-magnetic exchange interaction across the heterointerface 1 " 4 . However, ferro- 
electric BFO film, if not fully clamped by substrate, would generate voltage-induced strains that could be 
transferred to the magnetic thin film across the heterointerface and modulate magnetism together with exchange 
interaction through magnetoelastic coupling. Influence of this strain on voltage-modulated magnetism in BFO- 
based heterostructures has remained largely unexplored since raised by Mathur 5 , though a giant voltage-induced 
strain of over 5% has later been observed in BFO thin films with mixed rhombohedral and tetragonal phases 6 ' 7 . 

In this article, we explore how strain affects the voltage-modulated magnetism in BFO-based heterostructures 
by taking the Co 0 9 Fe 0 !(CoFe)/BFO thin-film heterostructure as an example. (001) BFO thin films were grown on 
(110) DySc0 3 substrate, exhibiting two-variant ferroelectric domains with 71° wall due to anisotropic film- 
substrate misfit strains 8 . In particular, the magnetic domain patterns in the top CoFe film almost copy the in- 
plane projection patterns of the ferroelectric domains at the BFO surface 2 ' 9 . Such magnetic domain pattern is 
induced by an effective magnetic field from the canted magnetic moment M c in BFO via Dzyaloshinskii-Moriya 
(DM) exchange interaction 1011 , namely, H DM -field. The direction of H DM -field (also M c ) is perpendicular to the 
plane of the polarization P and the antiferromagnetic axis L 1213 , i.e., H DM =P X L, resulting in a non-uniform 
H DM -field distribution based on the two-variant 71° ferroelectric domains. In this case, when electrically rearran- 
ging ferroelectric domain configuration to switch the in-plane net polarization in BFO by 180°, rotation of H DM - 
field within individual ferroelectric domain could lead to an overall 180° switching of the in-plane net magnet- 
ization 2 in CoFe. Detailed experimental analysis 14 reveals that the H DM -field (M c ) lies along the [1 12] and [112] 
directions corresponding to the [III] and [111] polarizations, while along the [ll2]and [112] directions corres- 
ponding to the [111] and [111] polarizations, respectively (figure la). 

Furthermore, the first few unit cells at the BFO thin film surface would become stress-free when the thickness of 
the BFO film was above a certain critical value [t c —70 nm for the BFO films grown on DySc0 3 (DSO) sub- 
strate 1516 ] to allow the presence of lattice defects such as dislocation (figure lb). In this case, sizable and non- 
uniform strain can be generated associated with local non-180 0 ferroelectric polarization (i.e., ferroelastic) 
switching under an electric-field. The strain can be further transferred to the top magnetic thin film across the 
interface and modulate the magnetic domain structure locally via magnetoelastic coupling 17 . Transfer of such 
non-uniform strain has recently been demonstrated in BaTi0 3 single crystal-based heterostructues 18 " 21 , which 
exhibit a similar one-to-one domain pattern match between magnetic thin film and ferroelectric BaTi0 3 under- 
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Figure 1 | Voltage-controlled magnetism in CoFe/BiFe0 3 heterostructure. (a) Locally coupled magnetic and ferroelectric domains in the 
Coo.9Fe 0 .i(CoFe)/BiFe0 3 (BFO) heterostructure reproduced by phase-field model (the first row), based on electric-field switching of the interfacial 
exchange coupling field H DM which is perpendicular to the plane of electric polarization P and the antiferromagnetic axis L in BFO (the second row). The 
slim arrows indicate orientations of local magnetization vectors, (b) Schematic of an epitaxial BFO film with partially relaxed substrate clamping when the 
film thickness t p exceeds the critical value t c for the generation of interfacial dislocations. 



neath. Particularly, it was proposed that 22 such non-uniform ferroe- 
lastic strain transfer alone could drive a 180° in-plane net magnet- 
ization reversal similarly to the non-uniform H DM -field driven 
reversal shown in figure la. These similarities thus raise one inter- 
esting question for BFO-based heterostructures: what role does the 
non-uniform strain in BFO play in the electric-field induced mag- 
netization reversal? A phase-field 23 ' 24 model is developed herein (see 
Method section) to address this question, and illustrate new possibil- 
ity in electric-field-controlled magnetism under coexistence of non- 
uniform H DM -field and non-uniform strain. 

Results 

In-plane net magnetization reversal purely by H DM -field. We first 
examine the influence of the interfacial H DM -field on the electric- 
field driven magnetization reversal in the CoFe/BFO heterostructure. 
If the electric-field applied along the [100] direction (£i 0 o) exceeds 
the coercive field of BFO (£ c ), the initial alternating [III] and [111] 
ferroelectric domains would switch by 71° to the [ill] and [111] 
domains, respectively, during which the polarization vectors would 
always keep a head to tail configuration (see figure 2a) to reduce the 
electrostatic energy. This leads to a full reversal of the average 
polarization along the [100] direction (i.e., Pi 0 o)> as illustrated by 
the simulated ferroelectric hysteresis loop in figure 2c. Note that 
individual ferroelectric domain under a certain £ 100 results in one 
unique local H DM -field distribution [Eq. (11)]. Correspondingly, 
figure 2b presents the in-plane projections of the H DM -field, e.g., 
an initial configuration of head to tail in-plane [TIO] and 
[TlO] orientations. The induced local magnetization distributions 
{viz. domain structures) in the CoFe film are almost identical to 
those of in-plane H DM -field (not shown here for simplicity), while 
both distributions are essentially the same as the in-plane projections 
of the ferroelectric domains and thereby accounts for the domain 
pattern transfer between the CoFe and BFO films 2 . Furthermore, a 
full reversal of the net [100] magnetization (M 100 ) in the CoFe film 
occurs when reversing the average in-plane H DM field along the [100] 
direction (H^) with an electric-field, as shown by their electric-field 
switching loops in figure 2e and 2d, respectively. The magnitude of 
the total H DM field [i.e., /z^ M , see Eqs. (11) and (12) in the Method 
section] is taken as 100 Oe based on a relevant experimental 
measumrent 9 . 

It is worth noting that such net magnetization reversal may not be 
triggered if the h^ M was too small, because the striped magnetic 
domains cannot be stabilized unless h^ M exceeds 35 Oe (see 



figure 3a), where the influence of strains is not incorporated 
(2 S =0). As it can be seen, once the striped domains form at 35 Oe, 
the net in-plane magnetization can then be reversed, i.e., the normal- 
ized magnetization M 10 o/M s changes from —0.7 to about 0.7, by 
changing the polarity of the net in-plane polarization. Simulated 
magnetic and ferroelectric domain structures at h^ M = 6 Oe and 
100 Oe are illustrated in the inset of figure 3a. We argue that the 
critical h^ M value for net magnetization reversal should be strongly 
dependent on the coercive field of the magnetic thin film, because 
influence of H DM -field is equivalent to an external non-uniform 
magnetic field [Eq. (9)]. 

Full magnetization reversal driven by both H DM -field and strain. 

Figure 3b shows the M 10 oAM s variations as a function of h^ M before 
(the upper panel) and after (the lower panel) reversing the net in- 
plane polarization, driven by both non-uniform H DM -field and non- 
uniform ferroelastic strain. Strain distributions before and after 
polarization reversal are plotted in figure 3c, with an alternating 
distribution of ±0.6% [also see Eq. (4a)]. Of interest, the magnetic 
striped domains can be stabilized during the growth of CoFe film by 
the growth-induced strains [see Eq. 4(b)] imposed by the striped 
ferroelectric (also ferroelastic) domain pattern at the BFO surface 
even when h^ M = 0 (i.e., no exchange interaction but sole strain 
effect, see figure 3b). However, these striped domains cannot be 
switched to achieve an in -plane net magnetization reversal when 
h^ M is small (weak exchange interaction). For example, at h^ M =6 
Oe, the striped domains only demonstrate a change in pattern 
periodicity (see corresponding domain structures in the inset of 
figure 3b) when reversing in-plane polarization. Thus, although 
non-uniform ferroelastic strain contributes to the stabilization of 
the domain stripes, it alone cannot trigger the net magnetization 
reversal. On the contrary, this strain could even suppress the 
reversal compared to the case purely driven by H DM -field 
(figure 3a), as demonstrated by (/) the enhancement of the critical 
value of h^ M from 35 Oe to 45 Oe (see the solid line in figure 3b), and 
{it) the different magnetic domain stripes at h^ M = 6 Oe (prevailing 
strain effect) and 1000 Oe (prevailing exchange interaction) with the 
same FE domain patterns underneath (i.e., both are before electric- 
field poling, see the upper panel of figure 3b). 

These counteractive non-uniform H DM -field and non-uniform 
strain can, however, be exploited to achieve a full rather than net 
magnetization reversal. As shown in figure 3b, the magnetic striped 
domains are straightened out to form a uniform single domain when 
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Figure 2 | In-plane net magnetization reversal via electric- field switchable H DM -field. In-plane electric- field (£ 100 /£ c = 1) induced dynamic changes in 
(a) the local distributions of polarization (top view) in the two-domain-variant BFO film, and (b) the in-plane projections of the H DM -field (the same for 
the local magnetization). The background color indicates the orientation of the local polarization/H DM /magnetization field (see the color wheel). 
Electric-field switching loops of average (c) polarization, (d) in-plane H DM field, and (e) magnetization along the [100] direction. E c represents the 
coercive field of the BFO film. 



h^ M is about 500 Oe. We attribute this to the almost equal contribu- 
tions from the competing strain effect and exchange interaction at 
this point. Indeed, taking the coupled magnetic and ferroelectric 
domain structures at h^ M = 0 Oe as the reference, the changes in 
magnetic free energy induced by strain and H DM -field are 8.26 X 
10 4 J/m 3 and -7.95 X 10 4 J/m 3 , respectively, at h° DM = 500 Oe. 
Furthermore, a full 180° magnetization reversal takes place when 
reversing the in-plane polarization (see corresponding domain struc- 
tures in figure 3b), which should promise higher application poten- 
tial for low-power memory applications 25,26 than the average 
magnetization reversal in striped multi-domain. Further simulations 
show that the range of h^ M required for such voltage-driven 180° 
magnetization reversal is from 393 Oe to 778 Oe, whereas the rest 
^dm vames mat are over 45 Oe only lead to a net magnetization 
reversal, in which the contribution form H DM -field is enough to 
trigger magnetization reversal but much smaller/larger than contri- 
bution from strain. Particularly, the range of h^ M for 180° reversal 
can be tuned by engineering the saturation magnetostriction coef- 
ficient ^ s (e.g., by using different magnets) and/or the magnitude of 
ferroelastic strains [e.g., by tuning polarization, see Eq. (4a)]. 
Detailed analysis can be found in Supplemental materials SI. 

Dynamics of magnetic domain switching. So far three different 
voltage-induced magnetic domain switching paths has been 
observed under the coexistence of H DM -field and strain, including 
the change in domain pattern periodicity (pattern exchange) when 



H DM -field is weak or absent (h^ M < 45 Oe), the full magnetization 
reversal when contributions from H DM -field and strain are 
comparable (393 Oe < h^ M < 778 Oe), and the net magnetization 
reversal for the remaining values of h^ M . To unravel the underlying 
physics, we track the time- dependent changes in the sum of the 
magnetostatic and exchange energy (i.e., F ms + F exch ) for these 
different paths (figure 4), and their local magnetic vector 
distributions at various time stages are shown in Supplemental 
Materials S2. For the net magnetization reversal at h^ M =1000 Oe, 
an energy barrier is surmounted during its initial stages (<0.04 ns), 
whereas for the change in domain pattern periodicity at h^ M = 6 Oe, 
a lower energy path is present throughout the evolution. Given the 
almost identical energy states between these striped domains with 
opposite net magnetization (i.e., pattern exchange and net reversal in 
figure 4), and also given their similar energy oscillation trends after 
surmounting the energy barrier, it can be concluded that the energy 
barrier can only be overcome by the unidirectional H DM -field rather 
than the uniaxial strain. In particular, for the full magnetization 
reversal at h^ M =500 Oe, the counteraction between H DM -field 
and strain leads to the lowest energy state before and after 
switching, however, an energy barrier still needs to be surmounted 
to complete the reversal. 

Discussion 

Voltage-controlled magnetism in magnetic/BiFe0 3 heterostructures 
is a complex process with multiple underlying mechanisms, depend- 
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Figure 3 | Effects of H DM -field and strain on voltage-induced magnetic 
domain switching. Variations on M 100 /M s (the normalized magnetization 
along the [100] direction) as a function of h^ M (the magnitude of the H DM - 
field) before and after reversing the in-plane average polarization from Pi 00 
to Py 00 in the CoFe/BFO heterostructure, driven by (a) only H DM -field, and 
(b) H DM -field together with ferroelastic strain. The insets are magnetic and 
ferroelectric domain structures at h^ M = 6 Oe and 100 Oe in (a), and the 6 
Oe, 500 Oe, and 1000 Oe in (b). The slim arrows indicate the orientations 
of local magnetization vectors, (c) Distribution of the non-uniform 
ferroelastic strain arising from the ferroelectric domains at the BFO surface 
before and after polarization reversal. 

ing on the size of such heterostructures 27 . In the present article, our 
focus has been put on the effects of both ferroelastic strains and DM 
exchange interaction in the CoFe/BiFe0 3 heterostructures currently 
investigated. Besides these effects, the modulation of magnetism 
could also be contributed by the voltage-induced changes in interface 
charge densities 28 ' 29 , and/or interface orbital configuration (e.g., via 
Fe-O hybridization, refs. 30,31). Such interface effects would become 
more remarkable as the thickness of the CoFe film (2.5 nm herein) 
further decreases 32 . On the other hand, further decrease in the lateral 
size of the heterostructures may allow us to study the complex inter- 
play among strain, exchange interaction, and possibly the charge/ 
orbital effect in ferroic single-domains. This is not only fundament- 
ally interesting, but also important for the design of high-density 
magnetoelectric devices 25 ' 26 . The model established herein could pro- 
vide a good starting point for the further study of these issues. 

In summary, a phase-field model has been developed to under- 
stand how strain affects the voltage-modulated magnetism in BFO- 
based heterostructures. The model is validated by reproducing the 
experimentally observed voltage- driven in-plane net magnetization 
reversal in a CoFe/BFO heterostructure 2 driven by non-uniform 
H DM -field. Then we evaluate the stable magnetic and ferroelectric 
domain structures by adding the influence of the non-uniform fer- 
roelastic strain that arises from the ferroelectric domain patterns at 
the BFO surface. Under such coexistence of H DM -field and strain, 
three different magnetic domain switching paths are discovered 
depending on the magnitude of the H DM -field, including a full 
180° reversal of uniform magnetization when contributions from 
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Figure 4 | Free energy analysis for different voltage-induced magnetic 
domain switching paths. Time-dependent changes in the sum of 
magnetostatic and exchange energy for three different voltage-induced 
magnetic domain switching in CoFe/BFO heterostructures, driven by both 
non-uniform H DM -field and non-uniform strain. The magnitudes of 
H DM -field (/2° M ) are taken as 6 Oe, 500 Oe, and 1000 Oe for pattern 
exchange, net magnetization reversal, and 180° full magnetization reversal, 
respectively. 

H DM -filed and strain are comparable. Analysis on magnetic domain 
switching dynamics demonstrates the lowest energy state for such 
full magnetization reversal, as well as the decisive role of H DM -field 
for both full and net magnetization reversals. 

Methods 

Phase-field model. In phase field approach, the magnetic and ferroelectric domain 
structures are described by the spatial distributions of the local magnetization vector 
M=M s (m 1 ,m 2 ,w 3 ) and local polarization vector P = (P 1 ,P 2 ,P 3 ), respectively, where 
M s and m f (i= 1,2,3) represent the saturation magnetization and the direction 
cosine 33 . 

The temporal evolution of the ferroelectric domain structure in (001) -oriented 
BFO thin films is governed by the time-dependent Landau- Ginzburg (TDGL) 
equation 34 , i.e., 



dPi(r,t) 
dt 



(1) 



where L is a kinetic coefficient that is related to the domain wall mobility and F P is the 
total free energy of the FE layer, which can be expressed as, 



(/bulk +/dastic +/electric +fpad)dV, 



(2) 



here/buik>/ekstic'/eiectric> and/ grad indicate the densities of bulk free energy, elastic 
energy, electrostatic energy, and gradient energy of the BFO, respectively, with V P 
representing the volume of the ferroelectric layer in the heterostructure. The math- 
ematical expressions for the/ b uik>/ g rad> and /electric of the BFO (001) thin films are 
described in literature 35 . Corresponding to the local in-plane electric-fields applied 
across the CoFe/BFO heterostructure via planar poling electrodes 2 , the electrostatic 
energy density /electric is obtained under a short-circuit boundary condition 36 . 



The elastic energy density / e ^ astic is calculated as, 



rP 
J elas 



P (r) 



(3) 



where c?. fc/ is the elastic stiffness tensor of the BFO and gy(r) the position-dependent 
elastic strain; gy(r) is the total strain and s y ~ p (r) is the spontaneous (stress-free) strain 
of the BFO arising from the electrostrictive effect, i.e., e y ~ p (r) = Q^/Pfc(r)P/(r) 
(i,j,k,l= 1,2,3,4), where Q y - fc / is the electrostrictive coefficient tensor. Following 
Khachaturyan's elastic theory 37 , the total strain e y (r) can be expressed as the sum of 
homogeneous and heterogeneous strains, i.e. e,y(r) = J^+d&ij(t). Among them, the 
heterogeneous strain d& v (r) does not cause any macroscopic deformation in a sample, 

i.e.,j<5£y(r)dV = 0, and can be calculated as 5eq(t) = [diii(r) / dxj + duj(r) / dxi] /% 

where «,-(r) is the displacement. 

The homogeneous strain 8[j represents the macroscopic deformation, whose in- 
plane components equal the film/substrate mismatch strain e\- if the BFO (001) thin 
films were fully constrained by the DSO (110) substrate, i.e., elT = e s n = —0.35%, 
&22 = e s 22 = —0.48%, and £12 =£2! = 0. The mechanical equilibrium equation, i.e., 
day/ dxi = ^(^/ekstic/^) / d x i = 0> * s men solved by taking ai 3 =a 23 = a 33 =0 38 . It is 
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noteworthy that such anisotropic biaxial in-plane strains can reduce the ferroelectric 
domain variants of the rhombohedral BFO from eight to two 8 , leading to the unique 
two-variant striped domains with 71° walls as reconstructed in figure la. Particularly, 
in a partly relaxed BFO thin film (figure lb), a non-uniform ferroelastic strain s^ s 
could be generated across the first few stress-free unit cells of individual ferroelectric 
domain, which can be expressed as, 



QmPkPidz. 



(4a) 



Here t s ( = t p -t c ) denotes the thickness of the stress-free unit cells at the BFO film 
surface (see figure lb), which is approximated to be 40 nm by taking the total 
thickness t p of BFO as 110 nm 39 . 

Corresponding to the two-variant FE domains at the BFO surface, a two-variant 
striped strain distribution with alternating ±0.6% can be derived. Such non-uniform 
strain e? s would further be transferred, at least partially 18 , to the upper CoFe film that 
has a much smaller thickness (t m —2.5 nm, ref. 2) than t s . Note that the two-variant 
striped domains in the BFO film should still emerge at the thickness t p of 1 10 nm as 
observed in experiments 39 , though the possible interfacial dislocations would some- 
what release the anisotropic mismatch strains from the orthorhombic DSO substrate. 
Nevertheless, even the fully constrained BFO thin films with striped domains can 
impose structural strains on the upper magnetic thin film during its growth, 33 simi- 
larly to those observed in magnetic/BaTi0 3 heterostructures. 18 ' 19 Thus, the average 
growth-induced strain affects the domain structure of the as-grown magnetic thin 
film together with the interfacial H DM field, which can be expressed as, 



growth _ 



- QmPtPjdz, 



(4b) 



where Pjj is the spontaneous (or remnant) polarizations under zero external electric 
fields. Accordingly, when applying an electric-field to the BFO film, the resultant 
polarization switching would change spatial distributions of both the H DM ( = P X L) 
field and the strain ejj s [Eq. (4a)] at the interface. 

Magnetic domain structures of the CoFe film will then evolve driven by these non- 
uniform H DM field and non-uniform £^ s (or the e| rowth for simulating the domain 
structures of as-grown CoFe film), and can be described by the Landau-Lifshitz- 
Gilbert (LLG) equation, i.e., 

'* ■ " ' ' M -y 0 MxH eff -— Mx(MxH eff ). (5) 



dz 



Here y 0 an d oc are the gyromagnetic ratio (taken as —2.2 X 10 5 m-A _1 -s _1 from ref. 
40) and the Gilbert damping constant (—0.01 from ref.41), respectively, whereby the 
real time step AT (—0.06 ps) for the magnetic domain evolution can be determined by 
Ar = i(l + a 2 ) /(y 0 Ms) with At = 0.02. H eff is the effective magnetic field, given as 
H e ff = — ( 1 / n 0 ) (5F m 1 8M) , with [i 0 denoting the vacuum permeability and F m the 
total free energy of the CoFe layer. The F m is formulated as, 



F M = | (fanis +/exch +/ms +/ H +/dastic) dV , 



(6) 



where / anis ,/exch>/ms>/H> and/ e ^ stic are the magnetocrystalline anisotropy energy 
density, exchange energy density, magnetostatic energy density, the H DM - field energy 
density, and elastic energy density, respectively. Among them, the/ anis is neglected for 
simplicity regarding the isotropic nature of the polycrystalline CoFe film. The iso- 
tropic /exch is determined by the gradient of local magnetization vectors, i.e., 



/exch =/ [(VmO 2 + (Vm 2 ) 2 + (Vm 3 ) 2 ] , 



(7) 



where / denotes the exchange stiffness constant. The magnetostatic energy density/ m; 
can be written as, 



fms 



-^Ms(Hf-m) 



yshape 



(8) 



(Hi). 



Here the stray field H^ ot consists of a heterogeneous part Hjj etero (m)from the mag- 
netostatic interaction which depends on local magnetization distributions and is 
obtained by solving magnetostatic equilibrium equation, 

i.e.,V • (^ 0 Hjj etero + ^ 0 M s m) =0, under periodic boundary conditions 42 . Hjf* also 
includes an demagnetization part H^ hape (m) that relates to the average magnetization 
m by the sample- size dependent demagnetizing factor matrix N 43 , i.e., 
Hf pe (m)=M s Nm. 

Furthermore, the H DM -field energy density can be expressed, similarly to the 
Zeeman energy of an external magnetic field, as, 



/ H =-^oMs(HS M .m). 
The H™ M indicates the H DM - field imposed on the CoFe film, and is given as, 



(9) 



,dz. 



(10) 



where t t denotes the thickness of the interface creating interfacial magnetic inter- 
action, and the H DM -field vector in the BFO layer can be expanded as, 

o I F\ Pi P3 \ I L\ L 2 L- 



h° 

= , = = ™ = = = (Pi,P2,Ps) x (P2P3, -P1P3, 0) (11) 
^P\+P 2 2 +Pl-^Pl{P\+P 2 2 ) K ) 

= , , ! *™, = , = (Pil* P 2 Pj, -P 2 P 3 -P 2 2 P 3 ). 

^p\+p 2 2+pi^pi{p\+piy 

The Pi and L t (i= 1,2,3) are the components of the polarization vectors and the 
antiferromagnetic axes along the cubic <100> axes, respectively, and the /z^ M 
represents the magnitude of the H DM field that depends on the chemical composition 
and geometric size of a specific BFO-based heterostructure 9 ' 44 . Note that the anti- 
ferromagnetic axis L is expressed as L=(Pi, P 2 , 0) X (0, 0, P 3 ), which indicates an 
alignment along either the in-plane [lTO] or [110] axis depending on polarization 
orientations (see figure la). Equation (11) clearly indicates that non-uniform nature 
of the H DM -field related to individual ferroelectric domain at the BFO surface, which 
can propagate across the heterointerface and act on the CoFe film similarly to the case 
in non-uniform ferroelastic strains sjj s [Eqs. (4a) and (4b)]. Combining Eq. (11), Eq. 
(10) can be rewritten as, 



h° 



DM f 
U J 



(P.P], P 2 P], -P\P3-P 2 2 P 3 ) 
V / P 2 +P 2 2 +P3 2 V^ 2 (^+^ 2 ) 



dz. 



(12) 



Influence of non-uniform ejj s on magnetic domain structures is characterized by 
changes in the elastic energy density of the CoFe, i.e.,/ e ^ stic , which can be calculated 
similarly to that in ferroelectric BFO [Eq. (3)] but use the elastic stiff constants of 
CoFe. Note these non-uniform sjj s are included in the position-dependent stress-free 
strain s°~ m (r) of the magnets, 



3/2^,^-1/3) (oref "*), (i=j) 
3/2A s m i m i + f , E |; s (or£? ow,h ), 



(13) 



with denoting the saturation magnetostrictive coefficient, taken as —62 ppm 
herein 45 . The coupling factor r\ (0 < rj < 1) is introduced to describe the possible loss 
of transferred strain due to imperfect interface contact, and is assumed to be 1 for a 
full strain transfer. 

Temporal evolutions of the ferroelectric and magnetic domain structures are 
obtained by numerically solving the TDGL and LLG equations using semi-implicit 
Fourier spectral method 46 and Gauss-Seidel projection method 47 , respectively. The 
material parameters used for simulations, including the Landau coefficients, elec- 
trostrictive coefficients, elastic constants, gradient energy coefficients of BFO layer, 
and the saturated magnetization, exchange stiffness constant, elastic constants of 
CoFe layer, are taken from literatures 48 " 50 . The discrete grid points of 128 Ax X 128Av 
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